Coupled Map Modeling for Cloud Dynamics 
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droplet, and advection. Through extensive simulations, the phases corresponding 



A coupled map model for cloud dynamics is proposed, which consists of the 



abatic expansion, fall of a droplet by gravity, descent flow dragged by the falling 



to stratus, cumulus, stratocumulus and cumulonimbus are found, with the change 

>< 

of the ground temperature and the moisture of the air. They are characterized by 
order parameters such as the cluster number, perimeter-to-area ratio of a cloud, and 
Kolmogorov-Sinai entropy. 

PACS# 47.52. +j, 92.40. Cy, 92.60. Jq 

Cloud dynamics plays an important role in the climate system, weather forecast, geo- 
physics and so on. However this elementary process in meteorology is very much compli- 
cated because it consists of different time and space scales and the phase transition from 
liquid to gas is coupled with the motion of atmosphere. Even if the flow of the atmo- 
sphere were known with accuracy using Navier-Stokes equation, we could not discuss the 
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morphology of cloud. In order to investigate such a complex system, construction of a phe- 
nomenological model is essential. In this letter, we introduce a coupled map lattice model 
of cloud formation which reproduces the diversity of cloud patterns. Characterizations of 
four cloud phases are also given. 

Coupled map lattices (CML) are useful to study the dynamics of spatially extended 
systems JTJ. Recently, CML has successfully been applied to spinodal decomposition 0, 
Rayleigh-Benard convection [EL the boiling transition M, and so on||. 

Here we construct a CML model of cloud in a 2-dimensional space. CML modeling is 
based on the separation and successive operation of procedures, which are represented as 
maps acting on a field variable on a lattice ||. Here, we choose a two dimensional square 
lattice (x, y) with y as a perpendicular direction, and assign the velocity field 
the mass of the vapor wl(x,y) and the liquid wj(x,y), and the internal energy E t {x,y) 
as field variables at time t. The dynamics of these field variables consists of Lagrangian 
and Eulerian parts. For the latter part, we adopt the following processes [J (1) heat 
diffusion (2) viscosity (3) buoyance force (4) the pressure term requiring divv to be 0, in an 
incompressible fluid []. Here, we use the discrete version of grad(divv) which refrains from 
the growth of divv. Indeed, we have already constructed the CML representations of the 
above four procedures [EJthat agree with experiments on the Rayleigh-Benard convection, 
which is a cardinal role for cloud dynamics. (5) diffusion of vapor wl(x,y); (6) adiabatic 
expansion; assuming the adiabatic process and the equilibrium ideal gas with gravity 
field, we adopt such an approximation that the temperature of the parcel risen from 
the height y to y + Ay is decreased in proportion to the displacement Ay. Thus the 
temperature E f (x,y) is decreased in proportion to Vy(x,y). (7) phase transition from 
vapor wl(x,y) to liquid w\(x, y) and vise versa accompanied by the latent heat. Here, 
we use the simplest type of the bulk water-continuity model in meteorology |J. The 
dynamics is represented as a relaxation to an equilibrium point w*(x, y) which is a function 

of temperature E t (x,y). (8) the dragging force; assuming that the droplets are uniform 

1 Here we consider such spatial scale that the centrifugal and Coriolis forces are neglected. 
2 If the vertical air motion is confined within a shallow layer, the motion of atmosphere can be regarded 
as incompressible flow |?J . 
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in size and fall with the terminal velocity V with neglect of the relaxation time to it, the 
dragging force is proportional to the product of the relative velocity (Vy(x,y) — V) and 
the density of droplets w\(x,y) at a lattice site. 

Combining these dynamics, the Eulerian part is written as the successive operations 
of the following mappings (hereafter we use the notation for discrete Laplacian operator: 
AA(x, y) = \{A(x -l,y) + A(x + l,y) + A(x, y - 1) + A(x, y + 1) - AA(x, y)} for any 
field variable A(x, y))\ For convenience, we represent state variables after an operation of 
each procedure with the superscript t + l/n where n is the total number of procedures. 

Buoyancy and dragging force 

< +1/3 (x, y) = vfa, y) + ^{E'ix + l, y) + E t {x- 1, y)-2E t {x, y)}-^w £ {x, y)(v y (x, y) - V) 
Viscosity and pressure effect 



V 



x, y) = if +1/3 (x, y) + vA^ 1 ^^, y) + ngrad(divv t+1/3 (x, y)) 



with grad(divv) as its discrete representation on the lattice ||. 
Thermal diffusion and adiabatic expansion 

E t+l l\x, y) = E\x, y) + \AE\x, y) - f3v l y (x, y) 

Diffusion of vapor 

wl +1/3 (x, y)) = w l v (x, y) + \Aw l v (x, y) 

Phase Transition To get the procedure, we use the discretized version of the following 
linear equations for the relaxation to equilibrium point w*\ 

- +a(w v (x, y) — w*) 



dw v (x, 


y) 


dt 




dwe(x, 


y) 


dt 




dE(x, 


y) 


dt 



a(w v (x,y) - w*) 

dw v (x,y) dw e (x,y). 



-Q{ 



dt dt 



Aexp(q/(E + const.)) if >W(x,y) 
W(x, y) otherwise 



which form is chosen to be consistent with the Clausius-Clapeyron's equation exp(- 
q/Temperature) |J, while W(x, y) = we(x, y) + w v (x, y) is the total mass of water. 

The Lagrangian scheme expresses the advection of velocity, temperature, liquid and 
vapor. This process is expressed by the motion of a quasi-particle on each lattice site 
(x, y) with velocity v(x,y). We adopt the method presented in ||, while for the liquid 
variable wg(x,y), we also include the fall of a droplet with a final speed V. Thus, the 
quasi-particle moves to (x+v x (x, y), y + v y (x, y) — V) to allocate wi(x, y) at its neighbors. 
Through this Lagrangian procedure, the energy and momentum are conserved. 

Summing up, our dynamics is given by successive applications of the following step; 
^(x.y) ] { ■u t+1 ' /3 (a;, y) ] [ -y t+2//3 (x, ] 

Tit , \ Buoyancy+Dragging I » I Viscosity+Pressure I j-, t , n I DifTusion 

< E\x,y) \ ^ 88 *l E\x,y) \ ^— ► I E\x,y) \ i — > 
. w\x,y) ) [ w\x,y) J [ w\x,y) J 

^ +2 / 3 (x,y) ) f ^+ 2 / 3 (x,i/) ) f ^ +1 (a;,y) 1 

J E t+1 / 3 (x,y) Ph -fc lsition EF-Wfay) Advecti ^ &a ^ ^(x.y) 
w* +1 / 3 (x, y) J [ w t+2 / 3 (x,y) J [ w t+1 (x,y) J 

For the boundary, we choose the following conditions; (1) Bottom plates: Assuming 
the correspondence between E and the temperature, we choose E(x,0) = E . (2) Top 
plates: We choose the no-flux condition E t (x, N y )—E t (x, N y — 1) = 0. For both the plates, 
we choose the no-slip condition for the velocity field, and adopt the reflection boundary 
for the Lagrangian procedure. The liquid and vapor are fixed at zero for both the plates. 
(3) Side walls at x = and x = N x : We use periodic boundary conditions. 

The basic parameters in our model are the temperature Eq at the ground, the Prandtl 
number (ratio of viscosity to heat diffusion z//A), adiabatic expansion rate (3, the terminal 
velocity of liquid droplets V, the coefficient for the dragging force 7, the phase transition 
rate a, the latent heat Q and the aspect ratio (N x /N y ). Hereafter we fix these parameters 
as A = 0.2, rj = v = 0.2, (3 = 0.2, V = 0.2, 7 = 0.2, a = 0.2, Q = 0.2, and study the change 
of the morphology in cloud as the ground temperature Eq and the total mass of water 
W = Y,x, y ( w v{.%iy) + w i{ x )V)) are varied (note that W is conserved). 

To see the spatiotemporal dynamics, the evolution of the mass of the liquid we{x,y) 
is studied. In Fig. [I], four typical time evolutions of Wi(x,y) are plotted. By changing E 
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Figure 1: 

Temporal evolution of cloud patterns. Snapshot of the mass of the liquid 
wj(x,y) is shown with the use of gray scale. The darker the pixel, the lower 
the liquid is. In other words, the white region corresponds to the high den- 
sity of liquids, i.e., to a cloud region. Snapshot patterns are plotted from 
top to down per 200 steps after initial 5000 steps of transients, starting from 
a random initial condition, (a) Stratus (E = 3.0, W = 0.006), (b) Cumu- 
lus (E = 3.0, W = 0.007), (c) Cumulonimbus (E = 4.0, W = 0.009), (d) 
Stratocumulus(£ = 5.0, W = 0.009). The lattice size is (N x , N y ) = 80 x 40. 



Stratus is a thin layered pattern of cloud, while cumulus is a thick lump of cloud. 
These two patterns are rather stable, while the other two patterns are dynamically unsta- 
ble. At stratocumulus, a thin layered cloud pattern is torn into pieces and small fragments 
of clouds are scattered. These scattered clouds vanish while a new layered cloud is formed 
again later. The formation and annihilation of clouds are periodically repeated. Cumu- 
lonimbus is a thicker cloud in height than cumulus. Besides the size change, the cloud 
pattern is unstable. The clouds split and coalesce repeatedly. The classification into four 
types is based on the comparison between our spatiotemporal pattern and the definition 



by meteorology ||, while the phase diagram is given in Fig. [2|, which is obtained from 
the pattern and quantifiers to be discussed. Summarizing the diagram, a cumulus or 
cumulonimbus is observed under the condition of rich moist air while a stratus appears 
in small W and under low temperature. 
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Figure 2: 

Phase diagram for the morphology of cloud. The term stratus, cumulus, cu- 
mulonimbus and stratocumulus correspond to the pattern (a),(b),(c) and (d) 
in Fig. [I], respectively. 



To classify these patterns quantitatively, we have measured several order parameters. 
First, we define a cloud cluster as connected lattice sites in which wg(x,y) is larger than 
a given threshold w c . f] C(t) is defined as the number of clusters disconnected with each 
other. Then, the "cloudiness" is measured by the total number of cloud sites, that is, 
S(t) = J2x=i^y=i®( w e( x iy) ~ w c)j where Q(x) is Heavisede function. The (temporal) 
average of the cluster number is large at the onset of cloud formation (i.e., small E and 
W), and at the stratocumulus. The quantity (S) / (C) measures the average size of each 

cloud cluster. It is larger at cumulus and is largest at cumulonimbus (Fig. 0-(a)). 

3 The phase diagram of the morphology of cloud does not depend on the choice of the threshold w c if 
0.02 <w c < 0.04. 
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Figure 3: 

(a) The average size of a cloud cluster (S) / (C) and (b) the stratus order 
parameter (SOP), plotted with the change of W and Eq. The average is taken 
over 20000 time steps after discarding 10000 steps of transients, starting from 
a random initial condition. 

To characterize the difference between stratus and cumulus, the morphology of clouds 
should be taken into account. Roughly speaking, the stratus is a one- dimensional like 
pattern while the cumulus is a two-dimensional one. To see the morphological difference, 
we have measured the total perimeter of cloud L(t) = Ylx=i!y=i J2s x =±i,6y=±i &{ w t{%, y) — 
w c )Q(w c —wi(x+5x,y+5y)). "Stratus order parameter" (SOP) is introduced as (l/(S(t)/(L(t)C(t)))) t , 
the inverse of the ratio of area to perimeter per cluster. If it is large the pattern is close 
to a one-dimensional object. Change of SOP with E and W is plotted in Fig. 0(b). As 
is expected by the definition of SOP and the thin nature of stratus cloud, it has a larger 
value at the stratus phase and takes a lower value at the cumulus. 

Of course, dynamical quantifiers are important to characterize the cloud patterns. For 
example, the fluctuations of the cluster size, SOP, and S/C are larger at the cumulonimbus 
and stratocumulus phases. To see the dynamics closely, we have also measured the time 
series of the spatial sum of the mass of the liquid L(t) = J2x, y w \{ x i y)-> which corresponds 
to the cloudiness of the total space. The evolution of Lit) is almost stationary at stratus 
and cumulus, with only tiny fluctuations. The change is periodic at stratocumulus and 
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chaotic at cumulus. It should be noted that the low- dimensional collective dynamics of 
the total liquid emerges even if the spatiotemporal dynamics is high-dimensional chaos. 

To characterize chaotic dynamics, Kolmogorov-Sinai(KS) entropy is estimated by the 
sum of positive Lyapunov exponents, as is plotted in Fig. f|(b). It has a larger value at 
stratocumulus and cumulonimbus, which implies that the cloud dynamics there is chaotic 
both spatially and temporally. It is also positive at a lower temperature that corresponds 
to the onset of cloud formation, where the dynamics is unstable. 




Figure 4: 

KS entropy calculated by the sum of positive Lyapunov exponents is plotted 
versus W and Eq. The first 20 Lyapunov exponents are computed by averaging 
over 20000 time steps after discarding initial 5000 steps. 

In summary, we have proposed a CML model for pattern formation of cloud by in- 
troducing a simple phase transition dynamics from liquid to vapor, so called the bulk 
water-continuity model |§. Our model reproduces the diversity of cloud patterns: stratus, 
stratocumulus, cumulus and cumulonimbus. This agreement implies that the qualitative 
feature of cloud dynamics is independent of microscopic details such as detailed droplet 
formation dynamics. 

In order to globally understand the phenomenology, the present computationally ef- 
ficient model is powerful, which makes us possible to characterize the cloud phases. We 
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believe that the observed phase diagram for the morphology of the cloud is valid for the 
cloud in nature. It is also interesting to propose that similar phase changes may be seen 
generally in convective dynamics including phase transitions, since detailed processes spe- 
cific only to clouds are abstracted in our model. Extensions of the present model to a 
three-dimensional case and inclusions of centrifugal and Coriolis forces are rather straight- 
forward. By these extensions, study of the global atmosphere dynamics of a planet will 
be possible. 

This work is partially supported by Grant-in-Aids for Scientific Research from the 
Ministry of Education, Science, and Culture of Japan, and by a cooperative research 
program at Institute for Statistical Mathematics. 
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